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As very-high-energy photons propagate through the extragalactic background light (EBL), they 
interact with the soft photons and initiate electromagnetic cascades of lower energy photons and 
electrons. The collective intensity of a cosmological population emitting at very-high energies (VHE) 
will be attenuated at the highest energies through interactions with the EBL and enhanced at lower 
energies by the resulting cascade. We calculate the cascade radiation created by VHE photons 
produced by blazars and investigate the effects of cascades on the collective intensity of blazars and 
the resulting effects on the extragalactic gamma-ray background. We find that cascade radiation 
greatly enhances the collective intensity from blazars at high energies before turning over due to 
attenuation. The prominence of the resulting features depends on the blazar gamma-ray luminosity 
function, spectral index distribution, and the model of the EBL. We additionally calculate the 
cascade radiation from the distinct spectral sub-populations of blazars, BL Lacertae objects (BL 
Lacs) and flat-spectrum radio quasars (FSRQs), finding that the collective intensity of BL Lacs 
is considerably more enhanced by cascade radiation than that of the FSRQs due to their harder 
spectra. As such, studies of the blazar contribution to the EGRB by Fermi will have profound 
implications for the nature of the EBL, the evolution of blazars, and blazar spectra. 



I. INTRODUCTION 

Since resolved blazars constitute the class of 
gamma-ray emitters with the largest number of iden- 
tified members 0, [|J , unresolved blazars are expected 
to have a significant contribution to the extragalac- 
tic gamma-ray background (EGRB) at GeV ener- 
gies [3|. The amount of this contribution remains 
unclear, ranging from a few to 100% of the EGRB 
pL S U 0, 11, S M HO, HI IH Q and depend- 
ing not only the distribution of blazars with respect 
to gamma-ray luminosity and redshift (the blazar 
gamma-ray luminosity function, GLF), but also on the 
distribution of blazar spectral indices (SID) at GeV 
energies [H, [IH, [HI, [l7| • Additionally, photons with en- 
ergies greater than 20 GeV can interact with the lower 
energy photons of the extragalactic background light 
(IR/Opt./UV; EBL) as they propagate from high red- 
shifts creating electromagnetic cascades of electrons, 
positrons, and photons. As blazars also constitute a 
cosmological population, the collective emission of un- 
resolved blazars at gamma-ray energies will exhibit an 
absorption feature at the high end of the Fermi LAT 
energy range. The shape of the absorption feature is 



sensitive to the blazar GLF and the EBL [lj|. The 
created cascade emission from absorption can enhance 
the collective blazar intensity at lower energies with 
the amount of enhancement being sensitive to blazar 
spectra at very-high energies [l9l |. 
II. CONTRIBUTION TO THE EGRB FROM 
CASCADES OF VHE PHOTONS FROM 
BLAZARS 

The contribution to the EGRB due to blazars can 
be viewed as the superposition of the collective inten- 
sity of intrinsic blazar spectra and the contribution 
from the cascade radiation from the interactions of 
very-high-energy photons from blazars with the EBL: 



I E (E ) = P E (E ) + r E (E ) , 



(f) 



where the intensity, Ie(Eq), is given in units of pho- 
tons per unit energy per unit time per unit area per 
unit solid angle emitted at observer-frame energy, Eq. 
The collective intensity of intrinsic blazar spectra in- 
cluding attenuation by the EBL is given by [for deriva- 
tion, seeLLSfl: 



Ih(Eo) 
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where Ef is some fiducial observer ray luminosity of a blazar at the fiducial rest frame 
frame energy (taken to be 100 MeV), energy, E'j — Ef(l + z'), p 7 = d 2 N / dL 7 dV corn is 
L 7 = EfdN(E' f )/dtdE is the isotropic gamma- the blazar GLF, t(E ,z) is the optical depth as 
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a function of observed photon energy and source 
redshift, Pi(a) = dN/da is the SID of gamma-ray 
blazars corrected for measurement error and sample 



bias [T|, 



1)(1 + Z) 01 EfFry 



-p7,min is the sensitivity of the survey above the 
fiducial observer frame energy 25]. In deriving the 
above equation, we have assumed the standard 
ACDM cosmology and that blazar energy spectra 
can be described as single power laws defined by the 
spectral index, a (dN 1 /dE 1 oc E~ a ), 

I 



Making the dependencies explicit, the cascade in- 
tensity is given by 



_ d 4 N c 
E ~ dtdAdVtdE 



(3) 



where dN c /dE is the spectrum of cascade photons due 
to pair production and Inverse Compton scattering 
[for derivation, seefl9j: 
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- T{E "' z) dE p dz, (4) 



where dNp/dE is the spectrum of Inverse Compton 
scattered radiation per electron of Lorentz factor, T, 
P(f;E p ,z) is the probability that the pair production 
interaction at a given redshift, z, of a primary photon 
of energy E p will produce electron-type particles of 
energies E e i = f x E p and E e2 = (1 - /) x E p , T 1 = 
f x E p /mc 2 , T 2 = (1 — /) xE p /mc 2 , and d 2 N 1 /dzdE p 
is the continuous spectrum of photons undergoing pair 
production interactions as a function of redshift and 
primary energy. Numerical integration of Eq. 0] is 
both complicated and computationally expensive, and 



thus, we make use of a Monte Carlo propagation code, 
Cascata which propagates photons from known 
spectra and redshift distributions and calculates the 
cascade spectrum of each photon. 

In determining the cascade contribution to the 
EGRB from VHE photons from blazars, primary pho- 
ton energies are randomly generated according to the 
absorbed EGRB intensity (in observer's quantities), 
which is the differential amount of radiation that has 
been "absorbed" by the soft photon background: 
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This intensity is interpreted as the radiation that cre- 
ates electromagnetic cascades. As such, it will also 
be used to normalize the cascade spectra in the final 
step. 



The absorbed EGRB intensity is also used to gen- 
erate the redshifts of the sources of the VHE photons 
by differentiating with respect to redshift: 



J 



dl% c 1 - e ~ T{ - E °' z * ) 

~dz~( Ea) = ~H~ 4nE 2 [n A + n m (l + z) 3 ]V2 J c 
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~E~ f 
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The created photons are then propagated until the 
redshift of interaction, randomly generated from the 
probability of interaction, p = 1 — exp[— t(Eq, z)]. 



III. RESULTS 



We determine the blazar contributions to the EGRB 
including cascade radiation assuming that the two 
blazar sub-populations, BL Lacertae-like objects (BL 
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Lacs) and flat spectrum radio quasars (FSRQs), form 
separate populations with respect to emission and evo- 
lution. We propagate photons from distinct GLFs and 
SIDs for the two blazar sub-populations. We assume 
BL Lacs to follow a luminosity-dependent density evo- 
lution model (LDDE) and FSRQs to follow the pure 
luminosity evolution model (PLE) [for more informa- 
tion about both models, see[l3j]. We determine the 
parameters of the SIDs for each blazar sub-population 
by applying the maximum-likelihood analysis of Ven- 
ters & Pavlidou [l6| to the spectral indices of the 
blazars listed in the Fermi LAT Bright AGN Sample 
Results for the Stecker et al. [H| EBL model 
and the Gilmore et al. 23] EBL model are plotted in 
Figures Q] and [2j respectively. 

In both EBL models, the effect of the cascade radi- 
ation is to flatten the overall collective intensity from 
blazars. Below 10 GeV, the cascade radiation for the 
two models is roughly comparable, though there is 
slightly more cascade radiation at the lowest energies 
for the Stecker model. Above 10 GeV, there is more 
cascade radiation for the Gilmore model than there 
is for the Stecker model due to the enhanced UV ab- 
sorption in the Stecker model. The resulting high- 
energy absorption feature in the overall spectrum is 
much more prominent in the Stecker case than in the 
Gilmore case owing to the higher UV background in 
the Stecker EBL model|26j. This is most especially 
true for the FSRQs, which follow the PLE model of 
evolution, resulting in more high-energy photons be- 
ing produced at higher redshifts than in the LDDE 
model [3. 

The hard spectra of BL Lacs cause the resulting 
cascade radiation to be considerable relative to that 
of the collective intensity from their intrinsic spectra 



even though their collective radiation tends to be con- 
centrated at low redshifts in the LDDE model[27[. In 
fact, the total intensity from BL Lacs could dominate 
the EGRB at high energies and even overproduce it, 
though it doesn't explain the lower energy intensity. 
By contrast, the relatively steep spectra of FSRQs are 
not conducive to much cascade radiation even though 
their collective radiation is concentrated more at high 
redshifts; hence, most of the collective intensity at- 
tributed to FSRQs occurs at lower energies. Thus, 
while the BL Lacs (including cascades) dominate the 
collective blazar intensity at the highest energies, the 
FSRQs dominate at the lowest energies. 

It should be noted that we have neglected the pos- 
sible effects of propagation through intervening mag- 
netic fields, which we will examine in a future publi- 
cation. 
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FIG. 1: Collective intensity of blazars including cascade radiation computed assuming separate SIDs and GLFs for BL 
Lacs (black; qo = 1.97, aa = 0.175; luminosity-dependent density evolution) and FSRQs (red; cto = 2.37, <ro = 0.14; 
pure luminosity evolution). Solid: Total collective intensity including both the collective intensity from intrinsic blazar 
spectra (including absorption) and cascade radiation. Dashed: Collective intensity assuming only absorption. Dot- 
dashed: Cascade intensity assuming the Stecker EBL model. Gray band: Preliminary Fermi EGRB (presented at the 
31 st International Cosmic Ray Conference by M. Ackermann). Data points: EGRET EGRB {triangles: Sreekumar et 
al. H; circles: Strong et al. HH) 
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surprisingly, models with high UV backgrounds will 
result in more suppression at high energies. 
[27] However, note that we did not take into account pos- 
sible breaks in the spectra. 
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FIG. 2: Same as in Figure [T] except for the Gilmore EBL model. 
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